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ABSTRACT 



43 
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Aims. With the afterglow of GRB 100621A being the brightest detected so far in X-rays, and superb GROND coverage in the 
optical/NIR during the first few hours, an observational verification of basic fireball predictions seemed possible. 
Methods. In order to constrain the broad-band spectral energy distribution of the afterglow of GRB 100621A, dedicated observations 
were performed in the optical/near-infrared with the 7-channel "Gamma-Ray Burst Optical and Near-infrared Detector" (GROND) at 
the 2.2m MPG/ESO telescope, in the sub-millimeter band with the large bolometer array LABOCA at APEX, and at radio frequencies 
with ATCA. Utilizing also Swift X-ray observations, we attempt an interpretation of the observational data within the fireball scenario. 
Results. The afterglow of GRB 100621A shows a very complex temporal as well as spectral evolution. We identify three different 

l/") ' emission components, the most spectacular one causing a sudden intensity jump about one hour after the prompt emission. The 

spectrum of this component is much steeper than the canonical afterglow. We interpret this component using the prescription of 

l/~) . Vlasis et al. (201 1) for a two-shell collision after the first shell has been decelerated by the circumburst medium. We use the fireball 

scenario to derive constraints on the microphysical parameters of the first shell. Long-term energy injection into a narrow jet seems to 
provide an adequate description. Another noteworthy result is the large (Ay = 3.6 mag) line-of-sight host extinction of the afterglow 
in an otherwise extremely blue host galaxy. 

Conclusions. Some GRB afterglows have shown complex features, and that of GRB 100621 A is another good example. Yet, detailed 
observational campaigns of the brightest afterglows promise to deepen our understanding of the formation of afterglows and the 
subsequent interaction with the circumburst medium. 



Key words, (stars) gamma-ray burst: general - (stars) gamma-ray burst: individual: GRB 100621A - Techniques: photometric 

1 . Introduction metallicity, dust) as well as to derive constraints on the progeni- 
tor (e.g., mass, rotation, binarity, supernova relation). 
1.1. The fireball scenario 

The late emission at X-ray to optical/radio wavelengths, the 
Gamma-Ray Bursts (GRBs) are generally accepted to be related S o-called afterglow, is dominated by synchrotron emission from 
to the death of massive stars. Due to their large gamma-ray lu- the external shock, i.e. e mission from relatiyistic electrons gy- 
minosity, GRBs can be detected to very high redshift, and thus rat i ng i n a magnetic field (Me szaros & Reesl[T997l:IWiiers et all 
provide a unique probe into the Early Universe. Understanding 1997; Wii ers & GaiamaL Il999h . This synchrotron shock model 
the emission mechanism and geometry is crucial for deriving the [ s widelyaccepted as the major radiation mechanism in the ex- 
burst energetics and number density, and observing and under- tern al shock, and the macroscopic properties of such shocks are 
standing the afterglow emission is of utmost importance to deci- largely understood. Under the implicit assumptions that the elec- 
pher the burst environmental properties (e.g., gas density profile, trons are "Fermi" accelerated at the relativistic shocks to a power 

law distribution with an index p upon acceleration, their dynam- 

* Based on data acquired with the Atacama Pathfinder Experiment ics can be expressed in terms of 4 main parameters: (1) the to- 

(APEX) under ESO programme 285.D-5035(A). tal internal energy in the shocked region as released in the ex- 
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plosion, (2) the electron density and radial profile of the sur- 
rounding medium, (3) fraction of the shock energy going into 
the ISM electrons e e , (4) the fraction of energy density in the 
magnetic field eg. Measuring the energetics and the energy par- 
tition (e e /eg) was possible only for a handful of the more than 
900 GRB afterglows so far, as it requires truly multi-wavelength 
observations between X-rays and radio frequencies. Moreover, 
there are large uncertainties in the microphysics: How are the 
relativistic particles accelerated? How is the magnetic field in 
the shocked region generated? What is its structure and evolu- 
tion? Addressing these questions is even more challenging. 

According to standard synchrotron theory, the radiation 
power of an electron with co-moving energy y e mc 2 is P e - 
4-/3cttc-)P;(B 2 /8tt), so that high energy electrons cool more 
rapidly. For a continuous injection of electrons, which is the case 
for ongoing plowing of the forward shock into the interstellar 
medium (ISM), there is a break in the electron spectrum, above 
which the spectrum is steepened due to cooling. This energy is 
time-dependent, so this frequency break moves to lower ener- 
gies for the ISM case and opposite for a wind medium. Since the 
spectral slope as well as the temporal decay slope are identical 
for the two density profiles, it is just the direction of the cool- 
ing break movement that allows to distinguish between ISM and 
wind density profile surrounding the GRB. 

Besides this cooling frequency v c , there is the injection fre- 
quency v,„, corresponding to the electrons accelerated in the 
shock to a power-law distribution with a minimum Lorentz 
factor, and the self-absorption frequency, v sa . The final GRB 
afterglow spectrum is thus a fou r -seg ment broken power law 
dMeszaros et all Il998t ISari et all 1 19981) separated by v sa , v„„ 
and v c . The order of v m and v e defines two types of spectra, 
namely the "slow cooling case" with v„, < v c , and the "fast 
cooling case" v m > v c . For each case , and depending on wind 
vs. ISM density profile, theory (Sari, 1999) predicts different 
slopes of the power law segments and speeds at which v m and 
v c should be moving. For "standard" parameters, v m should be 
moving from 10 14 Hz to 10 12 Hz within the first day, and v c 
from 10 17 Hz to 10 14 Hz. Due to sensitivity limitations in the 
sub-mm range, and lack of coordinated multi-wavelength obser- 
vations, there is not a single GRB data set sufficient (in terms 
of wavelength and temporal coverage) to unambiguously ver- 
ify these predictions for both frequencies, and just two GRBs 
where the high-frequency break (interpreted as cooling break) 
has been unambig uously shown to move dBlustin et all [2006; 
iFilgas et 



unambigt 
all 120121 . 



1.2. GROND and GRB 100621 A 

GR OND, a simultaneous 7-channel optical/near-infrared ima- 
ger (Greiner et al., 2008a) mounted at the 2.2 m telescope of the 
Max-Planck-Gesellschaft (MPG), operated by MPG and ESO 
(European Southern Observatory) at La Silla (Chile), started op- 
eration in May 2007. GROND has been built as a dedicated GRB 
follow-up instrument and has observed basically every GRB vis- 
ible from La Silla (weather allowing) since April 2008. The 
spectral energy distribution (SED) obtained with GROND be- 
tween 400-2400 nm allows us to not only fin d high-z candidates 
dGreiner et all l2009al iKriihler et all l2011al) but also mea sure 
the extinction and the power law slope dGreiner et all |201 II) . In 
the majority of all cases, this allows for a relatively accurate ex- 
trapolation of the SED into the sub-mm band, and consequently 
a prediction of the flux for sub-mm instruments, provided that v m 
has already passed the sub-mm band (which will be shown be- 
low to be the case for the majority of GRBs after about 1 day). 
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Fig. 1. GROND /'-band finding chart of GRB 100621 A, includ- 
ing the photometric comparison stars (roman and arabic letters). 
North is up, and East to the left. 



GRB 100621 A trig gered the Burst A lert Telescope (BAT) 
on the Swift satellite (Gehrels et al., 2004) on June 21, 2010 at 
T = 03:03:32 UT (Ukwattaet al., 2010a). The prompt emis- 
sion consists of a bright (25000 cts/s peak count rate in the 
15-350 keV band), smooth, triple-peak burst with a duration 
of nearly 70 s. Swift slewed immediately and started taking 
data with the XRT and UVOT telescopes at 76 s after the trig- 
ger. A bright X-ray afterglow was found at RA (2000.0) = 



2 1 h 1 m 1 3 S .24, Peel ( 2000.0) = -51°06'21'/7 with an error radius 
of l'/7 dEvans et allhoiOh . In fact, GRB 100621A has had the 
brightest X-ray afterglow ever detected: with an initial count rate 
in excess of ^140000 cts/s, it saturated the XRT CCD for sev- 
eral minutes. Starting 80 seconds after the burst, the X-ray light 
curve in the 0.3-10 keV band can be modelled with 4 power- 
rawfl. with decay indices and temporal breaks as follows: t?i = 
3.87+0.02, f breakl = 439+10 s, a 2 = 0.51!™ 2 , f break2 = 6.2^J 



ks, ff 3 = 1.0+-0.1, ?bre ak3 



122_ ; 2 i ks, 



and a 4 = 1.73+0.08 



(IUkwattaetalll2010bl) . 

GRB 100621A was also detected with INTEGRAL/SPI- 
ACSQ and Konus-Wind, providing a time-integrated spectrum 
with best-fit low-energy power law index -1. 7, high-energy in- 
dex - 2.45 and a peak energy £^=95* J 3 keV dGolenetskii et all 
2010). At z - 0.54 and standard cosmology (ff o =70 km/s/Mpc, 
Qm=0.27, Qa=0.73), this im plies an isotropic energy release of 
£ iso = (2.8+0.3) x 10 52 erg dGolenetskii et all 12010) . 

Initially, no UVOT counterpart was detected, and also rapid 
ground-based imag ing with robotic telescopes (like ROTSE, 
iPandev et all d20 101) 1 did not find an afterglow. Prompted 
by the discovery o f a very red afterglow with GROND 
(Upd ike et all 1201 Oi but see below), a spectrum taken with 
X-Shooter at the VLT determined a redshift of z=0.542 



1 Throughout this paper, we use the definition F„ oc t "v & where a is 
the temporal decay index, and fi is the spectral slope. 

2 http://www20 1 1 .mpe.mpg.de/gamma/instruments/integral/spi/acs/ 
grb/trigger/2010-06-21T03-03-26/index.html 
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cal/NIR imaging was calibrated against the primary SDSS0 stan- 
dard star network, or cataloged magnitudes of field stars from the 
SDSS in the case of g'r'i'z' observations or the 2MASS catalog 
for JHKs imaging. This results in typical absolute accuracies of 
±0.03 mag in g'r'i'z' and ±0.05 mag in JHKs ■ The light curve 
of the GRB 100621 A afterglow in all 7 GROND niters is shown 
in Fig. [2] 

2.2. Swift XRT data 

Swift/XRT data have been reduced using the XRT pipeline pro- 
vided by the Swift team. The X-ray spectra were flux-normalized 
to the epoch corresponding t o the GROND observat ions using 
the XRT light curves from lEvans et all (120071 120091) . We then 
combined XRT and Gala ctic foreground extinction (E(B -V) - 
0.03 mag: ISchlegel et alj d 19981) ) corrected GROND data to es- 
tablish broad-band spectral energy distributions (SEDs) which 
are shown in Fig. [3] 

2.3. NTT observations 

NTT/SOFI at La Silla was used to obtain NIR-spectroscopy. 
After recognizing the sharp drop in intensity at about T„ + 10 
ks we took four 60-s /-band images starting at 07:05 UT, on 
21 Jun 2010. While the results of the spectroscopy are deferred 
to a later publication (these are of no relevance for the purpose 
of this paper), the imaging provides an additional photometric 
data point at a time when no GROND observations were possi- 
ble anymore due to visitor mode regulations. The SOFI images 
were reduced in the same manner as the GROND JHK data (ac- 
tually within the same GROND pipeline), and calibrated against 
the 2MASS catalog. 



Fig. 2. Afterglow light curve of GRB 100621 A as observed with 
Swift in X-rays (top) and GROND in its seven filter bands (bot- 
tom). The /-band data points at 14 ks are from SOFI imaging, 
and the HK s -band data at 20 ks from a GROND-observation in 
morning twilight at which the /-band was already saturated by 
the rising Sun. The 7 vertical lines mark the times at which spec- 
tral energy distributions have been extracted (see text and Fig. [3}. 



( Milvang- Jensen et al., 2010), and also faint UVOT detections 
were recovered (Ukwatta et al., 2010b). 

Here, we describe our multi-wavelength observations and re- 
sults for GRB 100621 A, and present an analyses of the data in 
the framework of the fireball scenario. 



2. Observations 

2.1. GROND observations 

Some of the GROND data of this burst, in particular the /-band 
light curve and the host measurements, have already been re- 
ported in (Kri ihleret all l2011bl) . Here, we report the full data 
set, including the multi-band light curve, and the SED evolution. 
GROND exposures automatically started 230 s after the 
Swift trigger, one of the fastest reactions of GROND @2. 2 m so 
far. Simultaneous imaging in g'r'i'z' JHK S continued for 3.05 
hrs, and was resumed on nights 2, 4, and 10 after the burst. 
GROND data hav e been reduced in the standard manner using 
pyraf/IRAF ([Todvl[T99l lKtocu Yoldasetall [2008b). The opti- 



2.4. APEX observations 

Since the SED slope, even after extinction correction, was rather 
steep, the predicted sub-mm flux density of ^50 mJy at 1 day 
after the GRB led us to sub mit a DDT proposa l to ESO for ob- 
servations with LABOCA (Siring oet alll2009l) on the Atacama 
Pathfinder Experiment APEXO which was accepted at very short 
turn-around time. 

LABOCA, the "Large APEX Bolometer Camera", is an ar- 
ray of 295 composite bolometers. The system is optimized to 
work at the central frequency of 345 GHz with a bandwidth of 
about 60 GHz. 

The first APEX/LABOCA observation was obtained 1.08 
days after the GRB, leading to a clear detection. Two other addi- 
tional observations were performed at 2 days (another clear de- 
tection) and 4 days (upper limit only) after the GRB. This makes 
GRB 100621 A one of the rare cases with a sub-mm light curve 
(see section 5.3). All these observations were carried out in pho- 
tometry mode. 

Immediately after the first epoch observation (done in pho- 
tometry mode), we obtained at 5:32-6:26 UT a complementary 
observation of GRB 100621 A in mapping mode, for an expo- 
sure of 7x 420 s and reaching a lcr sensitivity of 14 mJy/beam. 
While no source was detected in this less sensitive observing 
mode, it verifies that there is no strong, unrelated source close to 



http://www.sdss.org 
4 APEX is a collaboration 



between the Max-Planck-Institut fur 
Radioastronomie, the European Southern Observatory and the Onsala 
Space Observatory. 
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Table 1. Secondary standards used for the GROND data 



Filter 



Star I 
21 01 12.58 
-51 05 17.2 



Star II 
2101 10.81 
-5104 54.6 



Star III 
21 01 15.88 
-51 06 17.4 



Star IV 
21 01 09.54 
-5106 22.2 



StarV 
21 01 05.82 
-5105 21.5 



Star VI 
21 01 08.30 
-5105 53.6 



Star VII 
21 01 14.38 
-51 05 25.8 



* 


16.60+0.05 


16.28±0.05 


18.54±0.05 


20.14±0.05 


20.34+0.06 


20.38+0.06 


19.49+0.05 


r' 


15.56±0.04 


15.64±0.04 


18.09±0.04 


18.58+0.04 


19.44+0.05 


19.70+0.04 


19.15+0.04 


i' 


15.29±0.04 


15.48±0.04 


18.00±0.04 


17.33±0.04 


19.18+0.04 


19.55+0.05 


19.10+0.04 


z' 


15.05±0.04 


15.31±0.04 


17.93±0.04 


16.69+0.04 


18.90+0.04 


19.39+0.04 


19.00+0.04 


J 


14.91±0.05 


15.34±0.05 


18.02±0.05 


16.23+0.05 


18.83+0.05 


19.42+0.05 


19.13+0.05 


H 


14.76±0.06 


15.31+0.06 


18.14±0.07 


16.04+0.06 


18.65+0.08 


19.42+0.09 


19.29+0.08 



Filter 



Star 1=1 

21 01 12.58 
-51 05 17.2 



Star 2 
2101 34.92 
-5105 59.3 



Star 3 
21 01 03.38 
-5103 26.6 



Star 4=IV 
21 01 09.53 
-5106 22.5 



Star 5 
210101.58 
-5107 43.8 



Star 6 
21 01 10.74 
-51 05 30.2 



K 



15.12+0.07 12.93+0.07 14.72+0.07 16.28+0.09 13.57+0.07 16.26+0.08 



Table 2. APEX/LABOCA observations at 345 GHz in photometry mode 



Date Time after GRB On+Off time Avg. t 
(UT) (days) (s) 


Flux EffNEFD 

(mJy ) (mJy sqrt(s)) 


Jun 22 04:38-05:30 1.0835 
Jun 23 07:27-08:15 2.1996 
Jun 25 07:51-08:42 4.2184 




607 
600 
592 


0.234 
0.358 
0.376 


35.5+3.3 61.8 

23.6+3.8 64.0 

5.2+3.4 54.4 


Table 3. ATCA observations 


Date Time after GRB 
(UT) (days) 


Flux @ 5.5 GHz 

G"Jy) 


Flux @ 9.0 GHz 

OjJy) 


Jun 24 19:00 -Jun 25 15:30 
Jun 25 15:30 -Jun 26 12:00 
Jul 17 08:00 -Jul 18 14:00 


4.0910 
4.9451 
26.2083 




137+17 
129+24 
-43+85 


150+28 
127+45 
49+100 



the GRB position, which otherwise could cause problems with 
the photometry mode data. 

Reduct ion of the ph otometric data was done with the soft- 
ware BoA (ISchullerL 120 1 2h using standard routines for photom- 
etry mode. Subscans were checked individually before averaging 
them together in order to identify and r emove outliers. The raster 
map was reduced with the CRUSH dKovacsL 2008) software 
package. Flux density calibration was done against Neptune, 
G45.1 andB13134. 



2.5. ATCA observations 



In response to the initial d etections of a bright afterglow 
of GRB 100621A dUkwatta et all l2010at lEvans et all [2010; 
lUpdike et all 120101: iMiTvang- Jensen et al., 2010), we also initi- 
ated observations of GRB 100621 A with the Australia Telescope 
Compact Array (ATCA) in Narrabri, Australia, at the frequen- 
cies of 5.5 and 9.0 GHz with an observing bandwidth of 2 GHz. 
The observation sessions were carried out between 24-26 June 
and 17-18 July 2010. The radio counterpart of the afterglow of 
GRB 100621 A was detected during the sessions carried out in 
June 2010 at both 5.5 and 9.0 GHz at a position coincident with 
those of the X-ray and optical counterparts, and it was unde- 
tected in the July 2010 session. 

It is possible that the observed decay between the first and 
second epoch, or part thereof, is due to interstellar scintillation, 
rather than due to the intrinsic decay of the afterglow. Otherwise, 
the fading at 5.5 GHz would have been rather early, indicating a 
rather low energy and/or eg. 



3. Overall light curve behaviour 



The overall temporal evolution of the afterglow at X-rays and 
the optical/NIR is shown in Fig. [2] The light curve in the X- 
ray band is very typical of X-ray afterglows as seen by Swift, 
with a steep decline (slope of -3...^-) during the first «400 s, 
followed by a shallow decay until about 122 ks, after which the 
decay steepens to a slope of 1.73+0.08 dUkwattaet al .. 2010b). 
In contrast, the temporal evolution of the optical/NIR afterglow 
is considerably more complex. From the start of the GROND 
exposures at 230 s post-trigger, the light curve shows a rapid 
rise with a\ - -4.3*i9. From about 400 s (consistent within 
errors with the end of the steep X-ray decline) to about 700 s, 
the light curve is more or less flat (a2 = 0.05 ± 0.05) with just 
a few wiggles. The sub-sequent decay has as - 1.15 + 0.15, 
significantly steeper than the X-ray decay at that time. After a 
short flattening (3^1 ks post-trigger), an extremely steep increase 
in optical/NIR brightness is observed from 4 to 5 ks after the 
trigger which has also been reported by the SIRIUS/IRSF team 
(Naito et al., 2010). This intensity jump is larger in the NIR than 
in the optical, reaching an amplitude of 1.9 mag in the /f,-band. 
A formal fit results in a\ - -14*H, the steepest flux rise we have 
ever seen in a GRB afterglow (at any time), both in the literature 
as well as in our GROND data over the last years. After a short- 
lived (5-9 ks) slow decline with a$ = 0.42 + 0.05, a steep decay 
with a6 = 2.3 ±0.1 sets in which flattens into the host flux level 
at around 3xl0 5 s. 
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4. Broad-band afterglow SED modelling 

4.1. Fitting framework, spectral breaks and cooling stage 

In the following, we will analyse our data in the framework of the 
fireball scenario, in particular in the formalism as described in 
iGranot & Sari (2002). From the single-epoch spectra in certain 
wavebands we can derive some basic boundary conditions, as 
follows. 

We start by fitting the GROND-data of the first 1 ks on its 
own. The SED built from the 7 GROND filters is very steep, but 
also clearly curved (right of center in Fig. |3), indicating substan- 
tial host-intrinsic extinction. As is standard practice, we apply a 
power law (as one segment of the fireball scenario) and fit the 
power law slope together with the dust extinction Ay in the rest- 
frame of the GRB (z=0.542). The resulting best-fit spectral slope 
in the optical/NIR range (well before the strong intensity jump 
at 4 ks) is measured to be yS ~ 0.8 + 0.1. Any slope flatter than 
ft ~ 0.7, in particu lar the theoret i cal p rediction of /3 - 0.5 for 
certain conditions dGranot & Sari, 2002), is safely excluded by 
the data (note that there is no ambiguity with the intrinsic host 
extinction Ay - 3.6 mag, see next section). 

Similarly, we fit the SwiftPiRT data on its own, and repro- 
du ce a slope of ft = 1.4 + 0.2 and Nh = 6.5 x 10 21 cm 2 as given 
in lUkwatta et al.1 (|2010b). Since we observe a steeper spectral 
slope in X-rays, this exclude s the fast cooling options ("spec- 
trum 4 & 5" in IGranot & Sari (2002)) at early times, and by 
construction (evolution from fast cooling to slow cooling) also 
at late times. 

Since the steepest possible fit to the GROND optical/NIR 
data is /? ~ 1.1, but the X-ray spectrum is significantly steeper 
than this, we are forced to introduce a break between the 
optical/NIR and X-ray data at intermediate times. Since at 
early times a single power law for the combined GROND and 
Swift/XRT data is sufficient, this break has moved into the cov- 
ered bandpass. We interprete this break as v c , as the observed 
slope difference of 0.6+0.2 is consistent with the predicted value 
of 0.5. If this break had moved from the infrared through the op- 
tical, the optical/NIR slope should have gotten bluer - which is 
not observed. In addition, the X-ray spectrum steepens, consis- 
tent with v c moving from high energies down through the X-ray 
band. We therefore conclude that the external density profile is 
constant (ISM-like). 

The simultaneous 5.5 and 9.0 GHz measurements at 4 and 5 
days after the GRB suggest a relatively flat slope of /3 ~ -0.25 
(with relatively larger error), implying that the self-absorption 
frequency v m is below 5.5 GHz. Again, as we observe (at cer- 
tain times) a spectral break between the optical/NIR and the 
X-ray bands, an interpretation according to "spectrum 2" or 
"spectrum 3" (Granot & Sari, 2002) with the self-absorption fre- 
quency slightly above 9.0 GHz (i.e. near its peak at the transition 
between y( 1 ~'')/ 2 /y~/ , / 2 to v 5 ^ 2 ) is impossible, as there would be no 
further break at higher frequencies. 

Thus, we are le ft with the option of "spectrum 1" 
(IGranot & SariL 120021) . for which the fireball prediction is f5 — 
-1/3 above the self-absorption frequency, in reasonable agree- 
ment with the measured/? = -0.25. While this conclusion is for- 
mally valid for the time of the radio measurements at 4 and 5 
days after the GRB , any other spectral ph ases ("spectrum 2" to 
"spectrum 5" from IGranot & Saril (12002b ') have been excluded 
by the above considerations. We therefore conclude that already 
at early times (T + 500 s) the afterglow is in the slow cooling 
phase. 

We therefore continue with the conceptual interpretation of 
slow cooling throughout our full data set, and the frequency or- 



dering as v sa < v m < v c , i.e. with the break between the op- 
tical/NIR and the X-ray part of the spectrum interpreted as the 
cooling break y c , and the break long-wards of the optical/NIR as 
the injection frequency v,„. 

We will model the SED at various epochs with a three- 
component power law, with slopes ft describing the radio range, 
ft the GROND range, and ft the X-ray range. According to the 
standard prescription (Granot & Sari, 2002), we fix the slope dif- 
ference to 0.5 around the cooling frequency v c , i.e. ft = ft + 0.5. 
We also fix ft = -1/3, due to the otherwise large effect on v m . The 
three power law segments are s moothly connected with a fixed 
smoothness parameter of 15 (see Beuerma nn et al.l ll999). 

4.2. Broad-band SED fitting 

For the following discussion, let us define 7 epochs which are 
sequentially in time: epoch 1 = 450-600 s (diagonal-hatched re- 
gion in Fig. |4]i epoch 2 = the sum of the time intervals 650-750 s, 
900-1 150, 1350-1800 s (cross-hatched regions in Fig.gji, epoch 
3 = 5.5-8.5 ks, epoch 4 = 94 ks, epoch 5 = 196 ks, epoch 6 
= 352 ks, epoch 7 = 416 ks, where the last three epochs are 
primarily determined by the times of the APEX and/or ATCA 
observations. In these latter three cases the optical flux has been 
determined by interpolating the GROND light curve which looks 
pretty smooth at these late times. The last three GROND epochs 
come with considerable systematic uncertainty due to the host 
subtraction. Due to the bright X-ray emission even at late times, 
no assumptions on the slope of the X-ray spectrum had to been 
made. 

A fit of these seven SEDs with the assumptions as listed 
at the end of the previous section and using all the available 
data at a given epoch is shown in Fig. [3] The most obvious re- 
sult is that the injection frequency (and there are good reasons 
why this is not a different break frequency, see above) moves 
to higher frequencies between epoch 5 (196 ks) and 6 (352 ks). 
This evolution is inconsistent with any prediction of the fireball 
scenario. While this is not a reason to condemn the fireball sce- 
nario, we discuss two possible options to explain this behaviour, 
both within the framework of the fireball scenario: 

(1) If one relaxes the usual assumption that the microphysi- 
cal parameters are constant, the break frequencie s would follow 
a mor e complicated evolution than described in IGranot & Saril 
(2002). While such a recourse has been offered for the descrip- 
tion of selected GRBs (e.g. Filgas et al., 2012), in the present 
case one would have to invoke an increase of e e proportional to 
r, or of eB as fast as f 3 ...f 4 . Moreover, this temporal evolution 
would be required only for the time between epochs 5 and 6, but 
not for the evolution as seen between epoch 4 to 5, or 6 to 7. 
Thus, we consider this option physically implausible. 

(2) Another option is that the true model, which results in the 
determination of the break frequencies, contains two (or more) 
different emission components which dominate at different fre- 
quency bands, or at different times. Already relatively small 
changes in flux of one component would lead to substantial 
changes in the break frequencies, even at constant slopes. A good 
example in our case is the epoch 3: assigning either all observed 
X-ray flux or just 50% of it (because the other 50% might be 
the normal underlying afterglow) to the component which pro- 
duces the large intensity jump in the optical/NIR will change the 
best-fit cooling break frequency by one order of magnitude. 

Thus, we conclude that a model-independent analysis of 
our data set is largely impossible, despite the broad frequency 
coverage and the multiple epochs available in all frequency 
bands. Moreover, as described above, the behaviour of the GRB 
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Fig. 3. Multi-epoch SEDs (different colours) of the late-time afterglow of GRB 100621 A as measured by Swift/XRT (right; not Nh- 
corrected), GROND (middle; not Ay-corrected)), APEX/LABOCA (middle left) and ATCA (far left), together with a broad-band 
model which fits all data available for the given epoch. The times of these SEDs are marked with vertical lines in Fig. and the 
resulting break energies given in Tab. |4] Since the optical/NIR and X-ray fluxes in epochs 1-3 are very similar, epoch 3 (jump 
component at 6.8 ks) has been scaled upwards by a factor of 20, and epoch 2 (flares) down-scaled by a factor of 4. The curvature 
in the GROND data is due to strong extinction of the afterglow light in the host galaxy (dotted line). The breaks seemingly show 
erratic variations in frequency - see text for an interpretation. Note in particular, that we consider the fits in this plot not to be the 
final physical interpretation of the data, as it links emission components at different wavelength regions which we argue in the text 
to not belong together. 



Table 4. Break energies as derived from the SED fitting as shown 
in Fig.[3j but see text. 



SED epoch 


Time 

(ks) 


(meV) 


Vc 

(keV) 


1 


0.52 


<550 


>8 


2 
3 


0.7* 
6.8 


<550 
<550 


9 n+0.6 
>8 


4 
5 


94 
196 


7 Q+0.9 

<2.0 


0.003-0.1 
0.003-0.1 


6 
T* 


352 
416 


10.6£j 


0.003-0.1 
0.003-0.1 



* This is the center of the first of three intervals - see text. 
** For these two epochs, our formal fit values for v„, are considered 
unphysical and thus likely an indication that the radio and optical/NIR 
emission stem from different components - see text. 



100621 A afterglow is so complex that we are also not able to 
test some predictions of (for example) the fireball scenario by 
our multi-epoch SEDs. 



Instead, the only approach left is to develop an interpretation 
as simple as possible within a given framework (and we chose 
the fireball scenario for this) which describes the data to a large 
(possibly full) extent. In what follows we use our data together 
with some basic arguments derived from the fireball scenario to 
disentangle the complex behaviour of the GRB 100621 A after- 
glow into several different components, the sum of which ex- 
plain the observed features. Our driving principle was to min- 
imize the number of assumptions, as well as emission compo- 
nents. This is likely not a unique description, and a more sophis- 
ticated interpretation is not excluded. 

We consider three different components: 

(1) a canonical underlaying afterglow, 

(2) flares during the first 1000 s, and 

(3) a jump component, most prominently visible in the opti- 
cal/NIR at 5.5-8.5 ks. 

Each of these components is allowed to have a different electron 
distribution p, and a different set of microphysical parameters 
such that the break frequencies in each are different. For most of 
the time, at least two of these three components overlap, and care 
has to be taken to assess which of the components dominates at 
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which time or in which spectral range. Our results, discussed be- 
low, suggest the following superposition of components, where 
the break frequencies are given for the dominant component in 
that frequency band: 

- epoch 1 : optical/NIR and X-rays dominated by canonical af- 
terglow, sub-mm and radio unconstrained; neither v c nor v,„ 
for SED of canonical afterglow are covered. 

- epoch 2: optical/NIR dominated by flares, X-rays are super- 
position with canonical afterglow, sub-mm and radio uncon- 
strained; 

- epoch 3: optical/NIR dominated by jump component, X-rays 
are ^50:50 superposition of canonical afterglow and jump 
component, sub-mm and radio unconstrained; neither v c nor 
v m of jump component covered. 

- epoch 4/5: optical/NIR dominated by jump component, X- 
rays dominated by canonical afterglow, sub-mm is likely the 
jump component; v„, of jump component in sub-mm. 

- epoch 6/7: optical/NIR still dominated by jump component, 
X-rays and radio dominated by canonical afterglow, sub-mm 
not constrained; v m of the SED of canonical afterglow is in 
the radio. 



4.2.1. Epochs 1 and 2 

At first glan ce, the rise time in the optica l is too fast for a for- 
ward shock dPanaitescu & Vestrandll2008l) . and the temporal and 
spectral parameters are not consistent with any closure relation 
(neither wind nor ISM density structure, with either standard or 
a jetted afterglow). Also, the subsequent part of the optical/NIR 
light curve (r o +300 to r o +600 s) is surprisingly flat. However, 
we note that the X-ray spectrum oscillates on a few hundred sec- 
onds timescale between a steep (J3 ~ 1.3) and a flat (J3 ~ 0.8) 
slope during the first few ks after the GRB. More interestingly, 
two of the three times of steep spectral slopes coincide with 
flux depressions in the (fluxed) X-ray light curve, and flux en- 
hancements (which could be described as optical flares) in the 
GROND data (lower panel in Fig.S at 300 and 700-800 s). This 
suggests that the evolution of the afterglow between r o +200 s to 
r o +2000 s is the superposition of two components, a "normal" 
afterglow and a flare component. 

In order to disentangle these two components, we fit the X- 
ray spectral index evolution (lower panel of Fig. |4} with a con- 
stant plus a number of separate Gaussians, whenever the spectral 
index deviates more than 3cr from the constant. We then apply 
a model composed of the rise and decline of a forward shock 
and the multiple Gaussians as derived in the previous step to the 
GROND light curve, now with fixed times of occurrence of the 
Gaussian components, but allowing different width and normal- 
izations. Due to better temporal resolution and S/N -ratio we con- 
centrate on the JHK S data. The residuals of such a fit without the 
Gaussians, i.e. the best-fit Gaussians to the GROND light curve 
on top of the forward shock are overplotted over the X-ray slope 
variation in the lower panel of Fig. |4] While there is no perfect 
agreement in all slope-oscillations, there is a surprisingly tight 
coincidence in the first two, at r o +300 s and r o +700 s. The re- 
sults of this exercise are: 

- the early rise in the GROND light curve is likely dominated 
by a flare, making the rising slope of the light curve partic- 
ularly steep; when including a flare at r o +300 s in the fit, 
the rise of the "normal" afterglow in the on-axis case is con- 
sistent with t 2 , suggestive of the canonical forward shock. 
This is additional evidence for a constant density profile, as 




700 



1000 



Time after GRB (sec) 

Fig. 4. Comparison of the fluxed X-ray light curve at 10 keV 
(top panel), the GROND J (yellow), H (blue), K s (green) bands 
(middle), and in the bottom panel the photon index of the X- 
ray spectrum (black, left y-axis scale) and the residuals of the 
model fit (see text) to the GROND JHK data (color as in the 
middle panel, right y-axis scale). The diagonal-hatched region 
denotes epoch 1, and the cross-hatched regions epoch 2. The 
dashed vertical lines mark the maxima in the photon spectral 
index F (r = f3 + 1) to guide the eyes. 



the rise in a wind profile woul d be much slower (f 05 to ? L0 ) 
dPanaitescu & Vestran"dll2008l) . 

- The relatively flat light curve during the interval at 300-800 
s after the GRB trigger is due to the contribution, and likely 
superposition, of flares. Once subtracted, the decay of the 
standard afterglow is flatter, namely a = 0.69 ± 0.03, where 
a systematic error of ±0.05 should be added due to the am- 
biguity of choice of the strength and width of the flares. 

- The optical/NIR emission during the intervals r o +450 s - 
r o +600 s (and T o +2500 s - r„+3000 s) are the only times 
when GROND sees "normal" afterglow at early times. This 
corresponds to our definition of epoch 1 . A combined fit of 
the GROND and Swift/XRT data results in a single power 
law of/3 = 0.81 ±0.02 with no spectral break being preferred 
over a fit with a break. Taking the corresponding Galactic 
contributions into account, the best-fit rest-frame dust extinc- 
tion and effective hydrogen absorption are Ay = 3.65 ± 0.06 
mag, and Nh = (1.8 ± 0.3) x 10 22 cirT 2 . The inferred slope 
above v c would be/^ = 1.31, with v c > 8 keV, and the cor- 
responding electron spectral index p — 2.62 ± 0.04. 

- The peak of the forward shock is at 380±30 s, corresponding 
to an initial Lo rentz factor of 71 ±3 ( according to the new 
prescription of Ghirland a et al.l (1201 2l) which returns values 
about a factor two lo wer than the previously used ones like 
l lMolinari et al. 2007h . 

- The emission during the flares is much steeper in X-rays, 
with best-fit spectral slopes in the 1.2-1.8 range. A com- 
bined GROND and SwiftptRT fit results in the need of a 
spectral break (at ~3 ke V), with low- and high-energy power 
law slopes of /? 2 = 0.86 ± 0.06 and /3 3 = 1.36 ± 0.06 (with 
fixed A/3 = 0.5). It is interesting to note that the spectrum al- 
ternates four times during the first 1000 s between this steep 
flare spectrum and the flatter "normal" spectrum. 



Considering these results for the 
with a Q = 0.69 ± 0.06, a x = 0.74 : 



'normal" afterglow, i.e. 
0.02 (note that we de- 
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viate from lUkwattaetail (1201 Obi) in that we fit the r o +700 s 
to Jo+lOO ks interval with one straight power law, but omit 
the higher-flux portion at T +6 ks, see below and Fig. |6}, and 
P2 = /?3 = 0.81 + 0.02 with inferred p = 2.62, we find consis- 
tency in the optical/NIR and X-ray decay slopes, but also note 
that this is much flatter than one would expect with the canon- 
ical closure relations for a standard afterglow with the given p 
in either wind (a = 1 .72) or ISM (a - 1 .22) environment. This 
suggests some form of energy injection. If the addition of energy 
is a power law in the observer time, £,(< t) oc f, then the flat- 
tening is b y Aa = e * 1.41(0.91) fo r a ISM (wind) density profile 
at v < v c (JPanaitescu et all 12006*) . Thus, with e = 0.35 - 1, de- 
pending on the circumburst density structure, consistency could 
be reached. As we will show below, our data are not compatible 
with a wind medium, so we adopt an energy injection according 
to Ei(< t) oc f - 35 until T + 4 ks. 
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Fig. 5. Early part of the GROND J-band light curve with a 
(slightly stretched in time) model of the two-s hell collisions 
overplotted (case 4, Fig. 7 in IVlasis et al. (2011)). Though this 
model was not aimed at reproducing the behaviour of the GRB 
100621 A afterglow, the similarity of the rise, structure at the 
peak and the decay slope is striking. The early part of the model 
should be ignored, as it depends on the relative timing of the for- 
ward shock of the first shell, the ISM density and initial Lorentz 
factor. 



4.2.2. Epoch 3 - the intensity jump 

While the short interval of the steep rise between 4.0-4.5 ks 
after the trigger is not covered by the Swift/XRT due to Earth 
limb constraints, the time of the first optical peak including 
the following slow decay phase until T +8 ks is covered with 
SwiftPiRT observations, but shows only a marginal X-ray flux 
enhancement, on the order of 50% relative to earlier and later 
times. This is in full agreement with the chromaticity seen within 
the GROND band (after host subtraction and extinction correc- 
tion), where the flux enhancement ranges between 200% (0.8 
mag) in the g'-band and 570% (1.9 mag) in the K s band, imply- 
ing a very red/soft spectral shape. A combined GROND/XRT 
spectral fit of the overlapping time interval 5.5-8.5 ks returns 
a single power law as best fit with a slope of /3 - 0.98 + 0.02 
when fitting all X-ray flux, or /3 = 1.0 + 0.03 when fitting 
just 50% of the X-ray flux (under the assumption that the other 



50% belongs to the "normal" afterglow). Two notes are in or- 
der: First, the SED can also be fit with a broken power law, with 
the break somewhere between the GROND and the Swift/XRT 
data. However, the improvement in reduced x 2 is only marginal, 
so we adopt the simpler model. Consequently, we assume v c > 8 
keV in the following. Second, the above decomposition assumed 
similar spectral slopes, which cannot be proven unambiguously. 
However, if the X-ray spectrum of the jump component would 
have been steeper by 0.5, with correspondingly v c being between 
the GROND and SwiftPiRT ranges, then one would not have ex- 
pected to see any X-ray flux increase at all. 



Time since GRB trigger [d] 
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Fig. 6. X-ray light curve of the GRB 100621 A afterglow with a 
broken power law fit, ignoring the enhanced emission at 5-8 ks 
which we assign to the jump component (see subsection l4.2.2b . 
The decay slope from around 1 ks nicely continues until 80 ks, 
when it steepens to a - 1.54 + 0.06. 



The overall shape of the rise, short shallow decay and subse- 
quent fast decay i s very similar t o the b ehaviour of the afterglow 
of GRB 081029 dNardiniet all 1201 ll) . where an analogous be- 
havior has been associated with the intrinsic properties of the 
GRB and not to changes in the intervening dust cont ent. In the 
meant ime, but independent of these observations, IVlasis et all 
(201 1) have presented numerical simulations of the collision of 
an ultra-relativistic shell in a constant density environment with 
the external forward shock, which produce similar flare light 
curves: Fig. [5] shows their case 4 model (with 2° half-opening 
angle; from their Fig. 7) overplotted over the GROND J-band 
light curve. In this scenario, the fast rise occurs when a second 
shell reaches the back of the first, self-similar Blandford-McKee 
shell. The steepness and amplitude of the rise depend on the half- 
opening angle of the jet, the Lorentz factor of the two colliding 
shells, and likely more parameters like the energy, the occur- 
rence time relative to the jet break, and eg. A parameter study 
much more extensive than that in Vlasis et al. (2011) is needed 
in order to be able to derive some of these parameters (or ranges 
thereof) for GRB 100621 A. However, a qualitative conclusion 
would likely be that GRB 100621 A has a large Lorentz factor 
or a small half-opening angle, or both. Among the sample of a 
handful of GRBs showing such features (Greiner, 2011), GRB 
100621A shows the steepest rise in time: a formal fit with T at 
the GRB trigger results in a^ e — 14 (which due to its late ap- 
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pearance is also insensitive on any possible change in T of the 
fit)! 

According tto IVlasis et al.l (|2011), the rather flat part after 
the jump is then due to the merging of the two shells, the heating 
of which compensates the fading flux from the forward shock 
of the first shell. After the jump, the light curves should follow 
the predicted slopes for the normal, single forward shock, but 
at a higher intensity level due to the additional energy injection 
by the colliding second shell. While this is difficult to convinc- 
ingly test with our data since the normal decay is not accurately 
enough constrained, the rise and the observed structure in the flat 
part of the light curve is surprisingly similar to the modelling in 
Vlasis et al. (2011), particularly their Fig. 6. We defer a more de- 
tailed comparison of this behaviour in GRB 100621 A with this 
shell-collision model to a future paper. 

4.2.3. The light curve beyond 20 ks 



We have shown in sub-section 14.2. 1 1 that the normal afterglow 
decay slope in the optical/NIR at a few ks after the GRB was 
ao - 0.69 + 0.06. An extrapolation of this decay at the same 
decay rate, i.e. with continued energy injection at the same tem- 
poral rate, underpredicts the later GROND data by at least a fac- 
tor of 2. Thus, the rate of energy injection would have had to 
increase over the early rate, if it were to explain the optical/NIR 
emission at r„+20 ks. We consider this unlikely, and thus con- 
clude that at late times, i.e. t > T o +20 ks, the optical/NIR fluxes 
are dominated by the process which led to the huge intensity 
jump at 4 ks. As the spectral shape of this emission was redder 
than that of the "canonical" afterglow, this statement will also 
be true for the sub-mm and radio bands (see next sub-section). 
At X-rays, we have shown in the previous sub-section that the 
contribution of the large intensity jump was marginal, at most 
50%, during the peak emission of the intensity jump. If the X- 
ray emission associated to the jump component subsequently 
dropped the same way as the optical emission, then it faded by a 
factor of 20 in the interval from r o +10 ks to T o +30 ks. The to- 
tal X-ray emission faded by just a factor of 2, implying that the 
X-ray emission beyond about r o +10 ks can be solely attributed 
to the normal afterglow. 

The fit to the X-ray light curve, using a broken power law and 
ignoring the enhanced emission at 5-8 ks, describes the overall 
behaviour very well. The break time is derived to be 80 ks, at 
which point the decay steepens to a — 1.54 + 0.06 (Fig.|6]l. 

This steepening of the light curve could be due to the cessa- 
tion of the energy injection. However, for our value of f3, a full 
cessation should lead to the canonical decay slopes of a — 1 .72 
(wind) or a — 1 .22 (ISM) in the standard afterglow scenario, or 
a - 1.96 or steeper for any jet model (see below). Thus, only a 
partial cessation of energy injection would be a viable solution. 

Alternatively, it could be the passage of the cooling break at 
continued energy injection. This would not work for the standard 
afterglow scenario of a spherical afterglow (i.e. F > 1/0, where 
6 is the jet half-opening angle), since the predicted slope change 
is just Aa - 0.25. Howev er, the predicted chang e is larger for 
a jetted outflow. Following Panaitescu et al. (2006), we consider 
two options: (i) a jet whose edge is visible and which does not 
expand laterally, and (ii) a jet with sharp edges which spreads 
lat erally and is observed w hen F x < 1. In their eqs. 34 and 
35. lPanaitescu et al. (2006) provide the flattening of light curves 
due to energy injection for the frequency range above and below 
v c . For option (i), the slopes depend on the circumburst medium 
density profile. Thus, we have three cases, each with a separate 
closure relation above and below v r . We start with the three cases 



for v < v c and determine e, the power of the energy injection (see 
above) such that the observed early decay slope of a - 0.72 is 
reproduced (we choose to take the value consistent with both our 
measured ao and ax, though this would not change our upcom- 
ing conclusion). With each of the three different values of e, we 
then check the predicted slope at v > v c for each of the three 
cases. Option (i) in the constant density environment returns the 
steepest slope, with a VK ^ = 1.2 (for e = 0.75). This is still flatter 
than the observed ax = 1.54 + 0.06 (Pig. |6j. The predicted slope 
depends only very weakly on f3, so also the trend of steepening 
fix towards the end of the observed X-ray light curve will not 
lead to consistency. We note that in this interpretation the energy 
injection is still active at the end of the X-ray light curve, i.e. at 
2 x 10 6 s, as we see no further steepening to a slope of a > 2.2 
(depending on any further softening of fix). 

Last, but not least, we note the coincidence of the measured 
slope of a — 1 .54+0.06 and the predicted ax — 1 .48 for the decay 
of the v > v c part of the afterglow in the spherical case. Thus, 
the steepening of the X-ray light curve at 80 ks could be due 
to the combination of both, cessation of energy injection AND 
passage of the cooling break in an ISM environment for an after- 
glow which is still in its spherical expansion phase (i.e. F > 1/0) 
when the collimation is not yet detectable. Admittedly, the need 
for this coincidence is not an attractive solution. At the moment, 
we have no more satisfactory explanation for the amount of the 
steepening of X-ray light curve at r o +80 ks. However, such 
break in the X-ray light curve is very common in the sample 
of ^700 Swift GRB afterglows, and thus a more generic problem 
dNousek et al.L 12006) rather than related to the specifics of GRB 
100621A. 



4.2.4. Epochs 4 and 5 

The two APEX/LABOCA detections correspond to a flux decay 



according to 



which then must accelerate considerably 



in order to be compatible with the upper limit at epoch 6. In 
the standard fireball scenario, the maximum in the sub-mm light 
curve is associated with the passage of the injection frequency. 
Since the observed decay slope is still considerably flatter than 
the expected t i(l ^ p ^ 4 for v < v m , the injection frequency of the 
dominating component should be near the LABOCA observing 
frequency during epochs 4 and 5. This is compatible with our 
best-fit SEDs: for epoch 4, the extrapolation of the GROND op- 
tical/NIR SED nearly exactly reproduces the APEX/LABOCA 
measurement, while for epoch 5 the optical/NIR flux (deter- 
mined from an interpolation between two GROND measure- 
ments) has faded more rapidly than the sub-mm flux, resulting in 
a move of v„, to lower frequencies. The speed of this frequency 
displacement between epoch 4 and 5 is measured as f- 115±0 - 55 ; 
consistent with the fireball prediction of f 3 ^ 2 . The observed op- 
tical/NIR flux is about a factor 2 above the extrapolation of the 
decay of the canonical afterglow, thus we assign this emission 
to the jump component. In contrast, as shown in the previous 
sub-section, the X-ray emission is due to the canonical afterglow 
component. Curiously, despite the steeper X-ray spectrum, a for- 
mal SED fit including the X-rays is possible due to the large gap 
between the optical and X-ray bands: since the X-ray spectrum 
has a steeper slope than the optical/NIR/sub-mm at this time, the 
large allowed range for v c can accommodate this bright X-ray 
component. 

Thus, with the two assumptions that (i) the contemporane- 
ously measured X-ray emission is a separate emission compo- 
nent, and therefore is left out from fitting; and (ii) the long- 
wavelength part is dominated by the jump component, we make 
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a combined spectral fit of epochs 3-6, where only epoch 3 con- 
tains X-ray data. We fix /?i = 1/3, A/3 - 0.5 between the 
GROND and the Swift/XRT band, and also fix the host extinction 
at the value of Ay = 3.65 mag as derived from the fit of epoch 3. 
With the lower S/N ratio of the later GROND SEDs, the slope in 
the GROND range is largely dominated by epoch 3, with a best- 
fit value for the combined fit of j3 2 = 0.90 + 0.04. The ATCA 
measurements then define the break energy v m as summarized 
in Tab. [4] We note that a fireball-compliant evolution of v„, from 
these values extrapolated backwards in time does not conflict the 

limit on v„, set by the NIR data at 5.5 ks (see dashed line labeled 
« r 3/2« in Fig rjj 

4.2.5. Epochs 6 and 7 

For these two epochs, we have the radio fluxes at two frequencies 
from the ATCA measurements. As described earlier, they are 
compatible with the v 1 ^ 3 slope as expected for the segment be- 
tween v ra and v m . At sub-mm, the APEX/LABOCA upper limit 
is well above this spectral component, and does not constrain the 
SED. 

The more or less unchanged radio flux in epochs 6 and 7 (for- 
mal fit results in t~ , though the large error bars of epoch 7 also 
allow a slightly rising flux) implies that the injection frequency 
is in the few GHz range (near our radio data). This conclusion 
is supported by two other observational constraints, namely that 
the radio spectral slope is somewhat flatter than v 1/3 , and that the 
radio flux must decline within the following 20 days in order to 
be compatible with the ATCA upper limits (see Tab. [3). 

A combined fit of the radio and optical/NIR data results in 
a best-fit injection frequency of order 2xl0 12 Hz, a factor 1000 
larger than our above estimate, and also a factor 10 larger than 
one would expect from the (fireball-compliant) evolution of the 
jump component. This suggests that the radio emission belongs 
to the canonical afterglow component, while the optical/NIR be- 
longs to the jump component (as argued above). This picture 
is consistent with a (again fireball-compliant) prediction of the 
early evolution of v,„, i.e. that v,„ is at frequencies shortward of 
the GROND NIR measurements at very early times (see the blue 
dashed line in Fig.fTOb. 

4.3. Characterization of the three emission components 

First, Tab. [5] summarizes the discussion from the above sub- 
sections with respect to the three emission components, rather 
than according to the epoch of observation, and Fig. |7]provides 
a visualization of the evolution of these three components with 
time. 



Table 5. Epochs at which the three emission components are 
seen at different frequency bands. 

canonical afterglow flares jump component 



Frequency (Hz) 



X-rays 


1, 3: partially, 4-7: fully 


2 


3: partially 


optical/NIR 


1: partially 


2 


3-7: fully 


sub-mm 


- 


- 


4+5: fully 


radio 


6+7: fully 


- 


- 




With these constraints on the varying combination of the 
three emission components at a given epoch, the combined fit- 
ting results in a total reduced xLa - El (162 for 145 degrees 
of freedom), thus being an acceptable fit. The best fit power law 



Observed energy (keV) 



Fig. 7. Visualization of the spectral energy distributions at the 
five epochs as discussed in the text, with panel 1 to 3 show- 
ing epochs 1-3, panel 4 showing epoch 4/5, and panel 5 showing 
epoch 6/7 . The frequency/energy ranges covered by our observa- 
tions are marked as shaded bands. Dashed lines mark the differ- 
ent emission components: afterglow (blue), flares (green), jump 
component (red). The thick line is the sum of these components. 



slope in the GROND range for the canonical afterglow (fully de- 
scribed in section 4.2.1) is /?2 = 0.82+0.02 with a strong host 
extinction of Ay = 3.65+0.06 mag, as already indicated by the 
very red colors of the afterglow. 

With the generic picture that the typical afterglow spectrum 
evolves from fast to slow cooling, we will now use the con- 
straints for each of the components, and try to infer a consis- 
tent picture of the evolution of the GRB 1 0062 1 A afterglow. The 
discussion is based on the formalism described in Granot & Saril 
(2002), and we use the same nomenclature of £52 = E/IQ 51 erg, 
andC= e e (p-2)l(p- 1). 



4.3.1. The canonical afterglow 

The SED of epoch 1 provides four constraints on the fireball 
parameters of the canonical afterglow: (i) a lower limit on the 
frequency of v c at that time (>8 keV), (ii) an upper limit on the 
flux density at v c (<0.035 mJy), (iii) an upper limit on v m based 
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on the non-detection of v,„ (or a break in general) in the GROND 
range, i.e. v m < 1.25 x 10 14 Hz (<2.4 /mi), and (iv) a lower 
limit on the flux at this limit frequency (> 9 mJy). Using the tw o 
equations each in lines 3 and 5 of Tab. 2 of Grano t & Saril (12002). 
these measurements translate into the following four conditions: 



(i) 
(ii) 
(iii) 
(iv) 



-3/2 

-1.62 



J2.12 



■ E'l 11 > 2.82 x 10 4 



,1.31 



£ e 



J/2 



E l 5 * 1 <7.17x 10~ 9 



Eli 2 < 6.46 x 10~ 6 



(1) 



el 1 -n 1/2 - £52 > 0.195 
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Fig. 8. Constraints on the microphysical parameters of the 
canonical afterglow component from the SED of epoch 1 at 520 
s after the GRB (black triangles), and of epoch 6 (colored ar- 
rows/lines). The limits on F(v m ) and v c from epoch 1 allow the 
parameter space above the lines of open triangles (left panel), 
and an upper limit of e e < 0.064 (top arrow on right panel). 
The limits from epoch 6/7 allow the parameter space below the 
lines of arrows (left panel), depending on the total energy. The 
requirement that e e < 1 translates into an upper limit in the den- 
sity and lower limit on e#, respectively (dotted line). The thick 
colored lines on the right panel show the corresponding allowed 
range for e e (where e e — e e x (p - 2)/(p - 1) = 0.39e e for the 
derived p=2.64). 

These equations define an upper limit one~< 0.064 (which 
translates into e e < 0.16 for our p-2.6A), and combined lower 
limits for eg and the external density as shown by the lines of 
arrows in Fig. [8] 

In principle, there are two more constrai nts, namely the lim it 
that the time of v„, crossing v c (5 — > 1 in dGranot & Saril 12002)) 
has occured within <520 s, and the transition 1 — > 2 (v m cross- 
ing v sa ) is constrained to >416 ks. However, these limits do not 
impose any additional constraints as shown in Fig. [8] 

Epoch 3 does not provide any further constraint on the 
canonical afterglow, as the X-ray flux and spectral shape can- 
not be independently differentiated from that of the jump com- 
ponent, as mentioned above. 

At epochs 4/5, the X-rays provide the only measurements of 
the canonical afterglow. Given the somewhat contrived conclu- 
sion that the late X-ray light curve after the break at 80 ks is 
due to a combination of cessation of energy injection and cool- 
ing break passage, and that it is still in the spherical phase, we 
refrain from adding these constraints here. 

Epochs 6/7, after re-fitting without the GROND optical/NIR 
data, provide no unambiguous measurement of v,„ and F(v m ) 
(or v c ), as the normalization of the power law segment which 



constrained. Stepping through v m in the range lxl0~ 7 keV to 
6xl0~ 5 keV reveals equally good fits as long as v,„ > 5 x 10~ 6 
keV (250 /mi). With this limit, we obtain: 



(i) 
(ii) 



e; 



4 /2 • Eg > 1.1 X 10- 3 



e' /2 -n 1/2 - £52 < 0.017 



(2) 



All combined constraints for the afterglow component are shown 
in Fig. |8] 



4.3.2. The early flares 

The flares are only observed at early times, and for a description 
we have picked epoch 2 to cover some of those. While we called 
these events flares, it seems obvious that these are somewhat dis- 
similar to the canonical X-ray flares observed by SwiftpCRT in a 
large fraction of GRBs: in the case of GRB 100621 A, the flares 
are prominent in the optical, rather than in X-ray s. If these have 
the sa me origin as the canonical X-ray flares (Margutti et al., 
2011), the only difference might be a lower peak energy. The 
broad-band spectrum between GROND and Swift/XRT is cer- 
tainly not a single power law (see section 4.2.1). As the low- 
energy part of a broken power law fit (j3 — 0.86) would be very 
steep for a Band function approach, the peak energy E pec ^ rather 
is below the GROND wavelengths, with possibly some exponen- 
tial cut-off at X-rays. Since the decomposition of normal after- 
glow component and flares is not unique, no statement can be 
made on a possible variation of the peak energy with time. 

4.3.3. The jump component 



10-' 
10-2 
10-3 
IO- 4 
10- 5 
10-6 




1 10 100 1000 10000 0.001 

ISM density n (cm- 3 ) 



0.01 0.1 



Fig. 9. Constraints on the microphysical parameters of the jump 
component, as derived from epoch 4/5 and 6/7 (see text) for the 
case of a constant ISM density profile. 



connects the 



,-1/3 



segment with the X-ray segment, is not 



As mentioned earlier, the IVlasis et all (1201 ll) interpretation 
of the optical/NIR emission at 5-8 ks is via the collision of two 
ultrarelativistic shells. 

The SED of the 5-8 ks event exhibits a straight power law 
of slope 1.1+0.1 covering the GROND optical/NIR and the 
Swift/XRT region. If interpreted using the lGranot & Saril (l2002h 
formalism for afterglows (the applicability of which is not ob- 
vious as the medium into which the colliding shell is evolving 
might be increasing in density, rather than being constant or de- 
creasing) the location of v c and v m remain ambiguous. If v c were 
longwards of 2.4 /jm (GROND Zf s -band), then the electron spec- 
tral index would be a reasonable p-2.2. However, in addition 
to the observed light curve decay at >10 ks being much steeper 
than the expected t~ , there would be further inconsistencies: 
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(1) if the circumburst environment had a wind density profile, v c 
would evolve to higher frequencies, i.e. into the GROND band, 
which is not observed; (2) if, alternatively, the circumburst den- 
sity profile were ISM-like, v c would move towards the LABOCA 
band. However, at epoch 4 the optical/NIR SED extrapolates 
nearly perfectly to the measured LABOCA flux, thereby not al- 
lowing any break. Thus, v c would have to be below 345 GHz at 
epoch 3. This would imply a later radio flux at least a factor 10 3 
larger than observed, and therefore can be excluded. We thus 
conclude that v c at epoch 3 must be >8 keV, implying a steep 
p=3.2. Using these constraints, we arrive at the following four 
conditions: 



(i) 
(ii) 
(iii) 

(iv) 



e -3/2 
-2.2 



■E: 1 ' 2 > 1.58 x 10 5 



-1/2 



J 52 
,1.6 



e 2 



.1/2 



E 2 $ l < 2.47 x 1(T 10 



E l J 2 < 2.36 x 1(T 4 



(3) 



ej 2 -n l/2 -E 52 > 0.358 



None of these conditions is violated by the constraints derived 
below for the emission of the 5-8 ks event. 

The SED of this 5-8 ks event is constrained by our measure- 
ments of epochs 4/5 and 6/7. During epoch 4, we measure v m 
and F(v m ), which provides the following two equations: 



(i) 
(ii) 



J/2 e-1/2 



4 /2 -« 1/2 



= 7.0x10" 



£52 =0.716 



"52 



(4) 



Epochs 6/7 provide an interesting constraint on the sub- 
mm/radio regime, despite the non-detections longward of the 
GROND-/: band. The APEX/LABOCA non-detection does not 
constrain the continuation of the optical/NIR slope into the nun- 
band, but a fireball-compliant extrapolation would suggest v m « 
1 x 10" Hz at epoch 6. Since we have argued earlier that the 
radio emission seen at this epoch at 5.5 and 9 GHz must belong 
to the canonical afterglow, we have to assume that the radio- 
component of the jump component must be self-absorbed to a 
level to not exceed the measured fluxes at 5.5 and 9 GHz. This 
results in v sa >0.8 x 10 11 Hz, i.e. v,„ = v s „ at epoch 6 (a nd 7) 
to within the errors. This is exactly what lVlasis et al.l(l201 lb find 
during the modelling of the radio light curve: the amplitude is 
strongly depressed due to self-absorption. 

In a constant external density profile, v sa is constant, and our 
above assumption does not violate any observational constraint 
at earlier or later times. For a wind environment, v sa decreases 
according to r 3//5 - this is slow enough that it does not conflict 
with the LABOCA detections at epochs 4/5. 

Thus, for the ISM case, we derive: 



(iii) 
(iv) 



•e 1/5 

2/5 



« 3 ^-£^ 5 =538 
7 ' 10 • E% m > 123 



(5) 



The combination of the last 4 equations translates into the two 
thin stripes of parameter space shown in Fig. [9] The resulting 
limits on the external density are rather high: since eg can- 
not be larger than 1, the external density must be ^20 cm" 3 . 
Moreover, the total energy is constrained to £52 > 0.2, and 
e e > 0.01. We stress again that these constraints are only valid 
if the lGranot & Saril (120021) formalism is applicable. 

5. Discussion 

5. 1 . Fitting assumptions and results 

The behaviour of the afterglow of GRB 100621 A at different 
epochs and frequencies has been found to be too complex rel- 
ative to our set of observational data to be able to constrain 



models. We therefore have adopted the fireball scenario and at- 
tempted to construct a consistent picture of the observed fea- 
tures. Before further discussion, we summarize our assumptions 
here: (i) First, we assume that the total emission is due to the su- 
perposition of 3 emission components; (ii) we have fixed A/3-0.5 
between X-ray and GROND power law slopes (whenever appli- 
cable); (iii) we have fixed p ra dw = -1/3 as derived from the two 
radio frequencies at epochs 6/7; and (iv) had to assume that the 
the jump component has to be self-absorbed in the radio. With 
these assumptions, we find a reasonably consistent picture which 
describes all of our observational facts (temporal and spectral 
slopes) except the slow X-ray decay at times >80 ks. 

For none of the three emission components in the afterglow 
of GRB 100621 A do we have enough observations at the right 
time to determine all fireball model parameters in a unique way. 
The constraints on these parameters as derived from our obser- 
vations are, in general, broadly consistent with expectations. The 
only inconsistent result is that for e^ of the afterglow component: 
the lower limits from epoch 1 are about 2 orders of magnitude 
higher than the upper limits as derived from epochs 6/7, assum- 
ing otherwise equal parameters (in particular total energy and 
density). There could be several reasons for this, one of which 
could be an evolving eg with time, though we do not consider 
this. A more obvious reason could be that the energy ejection 
(which was deduced to make spectral and temporal slopes in the 
early phases consistent with the fireball scenario) introduces a 
time-dependent variation between low- and high-frequency seg- 
ments (at radio wavelength, the impact of the energy injection 
will come later than at X-rays). This invalidates our assumption 
for epochs 6/7 in deriving constraints on v m , in that the radio 
and X-ray sections of the SED reflect the same internal energy 
budget. We therefore neglect the v,„ constraints from epochs 6/7 
in the following. If we allow v m to be just above 9 GHz during 
epochs 6/7, then no conflicting constraints are imposed anymore. 







1 


1 


1 1 -i- 1-+ _ ^ 1 1 


1 


I E 


10" 




X 


* 




— -ft— ■— _ 

a 


~-- --.^t- 1 /^ - — — . ! 


10 17 








G •-■ 


p 












co 


^ 


~~~~ — — - 1 










fi -..- 














CO 










10>« 








CJ 








! 


1015 




n 




^^ 








^ 



































« 










10" 









ft 


.1 


-s 












""-^ 


E 


in13 








""-v.^ 


" 


1U 1J 








*V>^_ 




'**^ 




= 


























'$0 "- 


























10 12 








"""•*. """-, 


"-x t-3/2 
















"■*, ■"">, 


• 






LABOCA 


"■"^ 




10" 









































103 



10 4 
Time since GRB (s) 



10 5 



106 



Fig. 10. Location of the two breaks v c (top end) and v,„ (bot- 
tom part) at different epochs in the late-time evolution of the 
afterglow of GRB 100621 A, for each of the three emission com- 
ponents (i) canonical afterglow (black), (ii) flares (red), and (iii) 
jump component (pink). Vertical bars indicate allowed ranges for 
v c or v,„. The wavelength coverage of our instruments is shown 
as vertical bars at the very left side. Dashed lines show the ex- 
pected evolution according to the standard fireball scenario af- 
ter obeying limits as derived from our observations at various 
epochs. 
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Despite the complex behaviour, we are able to unequivo- 
cally deduce a constant ISM-like circumburst density profile. 
The slow intensity decline of the external forward shock sug- 
gests continuous energy injection at a rate proportional to f 035 
during the first hour after the GRB. With the onset of the jump 
component, another sudden increase in energy happens which 
lifts the energy budget by a factor 2-5. 

One could imagine that the canonical afterglow and the early 
flares experience the same external ISM density, i.e. that they 
originate co-spatially. In this case, the combined constraints im- 
ply that the external density n >50 cirT 3 , otherwise the F{v m ) 
limit for the afterglow component would be violated. This in turn 
would imply that the energy driving the flares would be of order 
Ej SO (1 < £52 < 5), which is surprisingly large though not ex- 
ceptional. Correspondingly, we deduce 0.014 <e e < 0.064 and 
e B > 10- 4 . 

For the jump component, as mentioned above, we derived 
n >20 cm -3 , This is interesting as one could have imagined 
that this component originates in the wake of the afterglow, i.e. 
in a region cleared by the forward shock. Howev er, we cau- 
tion (again) that the interpretation with the Granot & Sari (2002) 
framework might not be appropriate at all. Further theoretical in- 
vestigation of such shell collisions are certainly warranted. 

5.2. Location of the dust 

From multiple SED fits during the early rise and early 
plateau (around 200^-00 s after the GRB trigger) we con- 
strain any variation of the extinction to AAy < 10%. The in- 
tense radiation of gamma-ray bursts has been repeatedly sug- 
gested to des troy the dust in i t s nea r environment through 
sublimation (Waxman & Draine, 2000; Fruch ter et all [2001; 
iPerna & LazzatiL 120021) . out to distances of a dozen parsec. The 
large dust column we observe in the afterglow of GRB 100621 A 
must therefore be at larger distances, most likely not related to 
the star formation site of the progenitor of GRB 100621 A. 
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Fig. 11. Comparison of our GRB 100621A sub-mm light curve 
to previous sub-mm observations of GRBs with more than 
one observation, and selected upper limits for a few famous 
GRBs. Different symbols mark different observer frequen- 
cies, and colors denote different GRBs ( except for the up- 
per limits) Data are from: GRB 030329: dKohno et all 12005b 
ISheth et all l2003h: GRB 090313: Greiner et al 20 13, in prep; 
GRB 080129: dGreiner et aill2009h : GRB 090423: dBocketall 
2009); GRBs 091102, 11 0709B , 110715A, 100901A, 110918A: 
ride Ugarte Postigo et afll2012l) . 



but for the majority ALMA will be the instrument of choice, 
once rapid turn-around target-of-opportunity observations will 
be offered. 



5.4. The GRB host 



5.3. Comparison with previous sub-mm detections 

Previous sub-mm meas urements of GRB a fterglows were 
initial ly non-detections dBremer et all [1998: Shephard et all 

998 J), and d etection s or even light curves are sparse 
dChandra et all 1200 8; Sheth et all 120031; iGreiner et all 120091; 

'erl ev et all 1201 2l IZauderer et all 120121) . Predictio ns of emis- 



sion at flux levels of several tens of mJy (e.g. Ilnoue et all 
2005) have not been materialized. So far, only a handful of 
GRBs have been detected in the mm/sub-mm, mostly using 
MAMBO at the IRAM 30m CChandr aeTatL 1200 8: Sheth et al., 
2003; Greiner etaill2009l). an dCARMA dChandra et al., 2007; 
Bock et al., 2009; Perle vet all 1201 2l) . GRB 100621Ais one of a 
handful of GRBs for which a sub-mm Tight curve' (more than 1 
detection) is available (Fig.fTTI). However, the complicated early 
optical/NIR light curve of GRB 100621 A makes even this rela- 
tively well-observed GRB too sparsely sampled in the sub-mm 
range, which leaves ambiguities in the interpretation of both, the 
light curve and the movement of the low-frequency break. 

Recent more aggressive attempts with APEX/LABOCA 
have confirmed to retu rn mostly non-detections 
ride Ugarte Postigo et all 120121) . indicating that the injec- 
tion frequency moves rather rapidly to frequencies below the 
LABOCA range, thus requiring sub-mm observations within the 
first day in order to achieve detections. APEX/LABOCA is able 
to do this for the best suited afterglows (steep optical/NIR SED), 



T he host galaxy of GR B 100621 A was extensively covered 
in iKriihler et all (1201 lbl) . including in addition to the GROND 
and Swift/UVCfT data. In short, the r'^21.5 mag galaxy is well 
detected from the UV (all Swift/UVOT filters) up to the /ri- 
band showing a very blue spectral energy distribution with 
(R - K)ab ~ 0.3 mag. The stellar population synthesis fitting 
of the host SED returns an age of the dominating stellar pop- 
ulation of only 0.05 Gyr, and an intrinsic extinction of A^? st = 
0.6+q .J ma 8' m stark contrast to the large afterglow (AG) extinc- 
tion of Ay G = 3.61+0.06 mag. The absolute magnitude of the 
host is Mb = -20.68 + 0.08 mag, and the star formation rate was 
determined as 13** M /yr. 

The APEX and ATCA non-detections of any flux at the po- 
sition of GRB 100621 A at >5 days after the GRB also provide 
first crude limits on the sub-mm and radio emission of the host 
galaxy, of <6.8 mJy at 345 GHz, <170/Jy at 5.5 GHz, and <200 
jt/Jy at 9 GHz (all 2<x confidence). Assuming that the dominant 
fraction of the radio emissio n would be of non -thermal origin, 
and using the formalism of lYun & Carillil (|2002), this implies an 
upper limit on the star formation rate of < 100 M /yr. 

Due to the bright, compact host, no observational attempt has 
been made with GROND to search for the supernova component 
which would have peaked about 6 magnitudes fainter (if extin- 
guished the same way as the afterglow) than the host brightness 
for a 1998bw-like SN-luminosity. 



13 



J. Greiner et al.: The unusual afterglow of the Gamma-Ray Burst 100621A 



6. Conclusions 

GRB 100621 A has shown the brightest X-ray emission after any 
gamma-ray bursts so far. Despite this, the afterglow at >200 s 
was not extraordinarily bright, and the strong host extinction 
made it only marginally detectable in Swift/UVOT observations. 
Yet, we obtained a decent data set with GROND as well as sup- 
porting APEX/LABOCA and ATCA measurements. 

The biggest surprise in the properties of the afterglow of 
GRB 100621A is undoubtly the sudden intensity jump after 
about 1 hr. Here, we have been able to characterize its proper- 
ties in hitherto unprecedented detail. The pecularity of this event 
is the complexity of the combined afterglow emission which we 
encounter. In order to disentangle this complexity, and to possi- 
bly even test afterglow models, a much denser sampling of the 
afterglow emission in time is required, both at sub-mm as well 
as radio frequencies. At least for sub-mm observations from the 
southern hemisphere, ALMA would be an ideal instrument if 
fast reaction times to external alerts like gamma-ray bursts can 
be implemented. 
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